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Abstract 

Several  previously  unreported  properties  of  turbulent  premixed  flames  were  measured  because  they  are  espe¬ 
cially  useful  for  the  future  assessment  of  direct  numerical  simulations  and  models.  These  new  properties  include 
local  stretch  rates,  a  wrinkling  parameter,  the  degree  of  flamelet  extinction,  and  the  reaction  layer  thickness,  which 
were  quantified  using  simultaneous  CH  planar  laser-induced  fluorescence/particle  image  velocimetry  (CH  PLIF- 
PIV)  diagnostics.  Other  reported  properties  that  are  useful  for  model  assessment  are  flame  surface  density  (Z)  and 
global  consumption  speed,  which  is  one  type  of  turbulent  burning  velocity.  Also  measured  was  the  Meneveau- 
Poinsot  stretch  efficiency  function  (/k),  which  plays  a  central  role  in  the  coherent  flamelet  model.  Some  images 
of  the  flame-eddy  interactions  show  how  eddies  exert  strain  and  how  flamelets  “merge.”  A  highly  wrinkled  (corru¬ 
gated)  flame  with  well-defined  boundary  conditions  was  stabilized  on  a  large  two-dimensional  slot  Bunsen  burner. 
It  was  found  that  the  turbulent  burning  velocity  of  Bunsen  flames  depends  on  the  mean  velocity  U,  which  was 
varied  independently  of  turbulence  intensity.  It  is  concluded  that  conventional  relations  for  the  turbulent  burning 
velocity  of  Bunsen  flames  are  inadequate  because  they  should  include  two  additional  parameters:  mean  velocity 
U  and  burner  width  W.  These  parameters  affect  the  residence  times  of  the  flame-eddy  interactions.  A  scaling 
analysis  is  presented  to  explain  the  observed  trends.  It  indicates  that  if  the  burner  width  is  sufficiently  large,  the 
long  flame  will  experience  significant  flamelet  merging,  which  is  one  factor  leading  to  the  “bending”  (nonlinear 
behavior)  of  the  burning  velocity  curve.  Images  of  CH  layers  show  that  flame  surface  area  is  lost  by  flamelet  merg¬ 
ing,  but  is  not  lost  due  to  local  extinction,  as  no  extinction  was  observed.  The  stretch  efficiency  function  increases 
with  increasing  integral  scale,  indicating  that  large  eddies  are  more  efficient  in  exerting  flame  stretch  than  small 
eddies. 
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Nomenclature 

Ap 

frontal  area  of  flame  brush  (Eq.  (5a)) 

U 

streamwise  velocity  averaged  over 

c 

mean  reactedness 

burner  exit  =  mR/ (pRLW) 

H 

height  of  Bunsen  flame 

U 

Favre- averaged  streamwise  velocity 

/o 

stretch  factor  (Eq.  (1)) 

UCL,0 

mean  velocity  at  origin  (Fig.  1) 

K 

stretch  rate  (Eq.  (9)) 

uf 

root-mean-squared  velocity  fluctuations 

Ka 

Karlovitz  number 

W 

burner  width  (Fig.  1) 

L 

length  of  burner  (Fig.  1) 

X 

streamwise  coordinate  (Fig.  1) 

mR 

mass  flow  rate  of  reactants 

y 

lateral  coordinate  (Fig.  1) 

M 

flamelet  merging  term  (Eq.  (2)) 

a  0 

thermal  diffusivity  of  reactants  at  296  K, 

m 

M/E 

0.18  cm2 

Ma 

Markstein  number 

<5ch 

CH  reaction  layer  thickness 

P 

perimeter  of  time-averaged  flame  posi¬ 
tion  in  laser  sheet  (Fig.  1) 

<$L0 

unstretched  laminar  flame  thickness, 

Ft 

perimeter  of  wrinkled  flame  surface,  av¬ 

=  7.4ofLo/^LO  —  0.35  mis 

eraged  over  70  images 

<$T 

turbulent  brush  thickness 

Q 

flamelet  quenching  term  (Eq.  (2)) 

0 

coordinate  normal  to  flame  brush 

SL0 

unstretched  laminar  burning  velocity  = 

stretch  efficiency  function  (Eq.  (3)) 

0.38  mis 

Vj 

turbulent  diffusivity 

S_T 

local  consumption  speed 

i 

integral  scale  of  turbulence 

Sj 

turbulent  burning  velocity  defined  as  a 

E 

flame  surface  density 

global  consumption  speed  (Eq.  (5c)) 

Q 

wrinkling  factor  (Eq.  (22)) 

1.  Introduction 

One  goal  of  this  work  was  to  quantify  several 
properties  of  turbulent  premixed  flames  that  have 
not  been  quantified  before  and  are  especially  use¬ 
ful  for  the  assessment  of  direct  numerical  simula¬ 
tions  (DNS)  and  models.  It  was  decided  to  select 
a  flame  that  is  in  the  “corrugated  flamelet  regime” 
and  use  simultaneous  CH  planar  laser-induced  fluo¬ 
rescence/particle  image  velocimetry  (CH  PLIF/PIY) 
diagnostics,  which  previously  had  been  applied  only 
to  nonpremixed  combustion.  These  diagnostics  al¬ 
low  one  to  identify  the  thin  wrinkled  reaction  zone 
(where  CH  exists)  and  the  nearby  velocity  field.  Prop¬ 
erties  that  were  measured  include  local  stretch  rates, 
a  wrinkling  parameter,  the  degree  of  flamelet  extinc¬ 
tion,  reaction  layer  thicknesses,  flame  surface  density 
(i7),  and  global  consumption  speed,  which  is  one 
type  of  turbulent  burning  velocity  that  is  described  be¬ 
low. 

Recently,  Bell  et  al.  [1]  obtained  impressive  DNS 
results  for  a  turbulent  premixed  flame  in  three  dimen¬ 
sions  with  complex  chemistry  and  an  adaptive  grid 
resolution  of  62.5  pm.  However,  to  determine  if  the 
grid  resolution  is  adequate  and  if  the  boundary  con¬ 
ditions  were  handled  properly,  comparisons  of  DNS 
results  with  measurements  are  required.  If  a  modeling 
approach  is  employed  (instead  of  DNS),  two  promis¬ 
ing  models  are  the  coherent  flamelet  model  (CFM), 
which  is  based  on  the  Reynolds-averaged  Navier- 


Stokes  and  flame  surface  density  equations,  and  the 
large  eddy  simulation  with  a  flame  surface  den¬ 
sity  subgrid  model  (LES-FSD).  The  CFM  has  been 
used  to  simulate  Bunsen  flames,  V-flames,  spherical 
flames,  and  planar  flames  by  Prasad  and  Gore  [2], 
Veynante  et  al.  [3],  Baritaud  et  al.  [4],  and  Duclos 
et  al.  [5],  respectively.  Alternatively,  the  LES-FSD 
model  [6-9]  has  been  applied  to  bluff-body  flames  by 
Knikker  et  al.  [8]  and  Weller  et  al.  [9]. 

To  assess  DNS,  LES,  or  CFM  results,  it  is  impor¬ 
tant  that  the  following  three  physical  mechanisms  are 
correctly  simulated. 

(1)  The  turbulent  burning  velocity  (Sr)  is  enhanced 
by  the  increase  in  surface  area  of  the  wrinkled  re¬ 
action  zone,  which  can  be  quantified  by  the  flame 
surface  density  (i7). 

(2)  This  increase  in  surface  area  of  the  reaction  zone 
is  caused  by  the  stretch  rate  ( K )  that  is  induced 
by  local  velocity  gradients  and  flame  curvature. 

(3)  The  stretch  rate  (K)  of  the  flame  is  caused  by 
vortical  motions  associated  with  turbulent  ed¬ 
dies. 

To  determine  if  these  mechanisms  are  simulated  cor¬ 
rectly,  measurements  of  the  quantities  Sj,  K,  and  E 
are  of  central  importance.  The  relationship  between 
these  three  quantities  is  shown  by  equations  that  were 
derived  by  Bray  [10],  Candel  and  Poinsot  [11],  and 
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Meneveau  and  Poinsot  [12]: 


—  ^LO^O  J  Zdr\, 

—oo 

~dz  ~ds  d2s 

U—  +  V—  =  vt— — ^  +  KS-Q-M, 
ox  oy  dyz 


K  =  ( u'/ODt . 


(1) 

(2) 

(3) 


Eq.  (1)  summarizes  physical  mechanism  (1):  the  tur¬ 
bulent  burning  velocity  is  proportional  to  the  wrinkled 
flame  area,  which  is  an  idea  that  was  proposed  by 
Damkohler  [13,14].  In  Eq.  (1)  the  flame  surface  den¬ 
sity  D  is  the  average  amount  of  flame  surface  area  per 
unit  volume,  and  r]  is  the  coordinate  that  is  normal  to 
the  flame  brush.  The  integral  of  E  drj  has  been  shown 
[1]  to  equal  the  ratio  Ax/Ap,  where  Aj  is  the  wrin¬ 
kled  flame  area  and  Ap  is  the  flame  area  if  there  was 
no  wrinkling,  /q  is  a  stretch  factor,  which  has  been 
shown  [15]  to  be  1-0.28  MaKa;  thus,  /q  is  approx¬ 
imately  1.0  if  the  fuel-air  ratio  is  chosen  to  yield  a 
Markstein  number  that  is  sufficiently  small,  as  was 
done  in  the  present  work.  Eq.  (2)  summarizes  physi¬ 
cal  mechanism  (2):  the  increase  in  flame  area  is  due 
to  the  mean  stretch  rate  ( K ),  as  well  as  the  quanti¬ 
ties  Q  and  M,  which  quantify  the  flamelet  quenching 
and  merging  rates,  respectively.  Eq.  (3)  defines  the 
Meneveau-Poinsot  stretch  efficiency  function  7~k  that 
relates  the  mean  stretch  rate  K  to  the  rms  velocity 
fluctuations  ur  and  the  integral  scale  i  [12].  Prior  to 
the  present  study  /k  could  be  estimated  only  from 
computational  and  experimental  studies  of  a  laminar 
flame  interacting  with  a  single  vortex.  The  preferred 
way  to  determine  /k  is  to  measure  the  stretch  rate 
within  turbulent  flames,  using  the  methods  described 
below. 


truly  universal  and  category-independent.  Therefore 
it  is  recommended  that  in  any  attempt  to  collapse  val¬ 
ues  of  turbulent  burning  velocity  to  a  set  of  curves, 
only  data  from  one  category  be  considered.  This  im¬ 
plies  that  there  will  be  a  different  plot  of  the  burning 
velocity  correlation  for  each  category. 

For  example,  two  parameters  that  control  Bunsen 
flames  are  the  mean  velocity  and  burner  width,  both 
of  which  affect  the  height  of  the  flame.  Taller  flames 
become  more  wrinkled  than  short  flames  because  it  is 
observed  that  wrinkling  increases  with  distance  from 
the  attachment  point.  Thus  it  can  be  expected  that  for 
Bunsen  flames,  a  burning  velocity  correlation  would 
depend  on 


ST  (  u'  l  U  W\ 

- =fcn[  - , - ,  Ma, - , -  I 

ShO  \SU)  <$L0  SU)  <$L0 / 


(4) 


where  Ma  is  the  Markstein  number  for  given  fuel  type 
and  equivalence  ratio.  Ma  has  been  shown  to  affect 
burning  velocity  even  for  very  high  Reynolds  condi¬ 
tions  [18].  The  unstretched  laminar  flame  thickness 
(<$Lo)  is  defined  as  the  temperature  difference  across 
a  laminar  flame  divided  by  the  maximum  tempera¬ 
ture  gradient;  this  distance  is  7.4  (apo/Spo)  based  on 
the  measurements  of  Bechtel  et  al.  [19].  For  spherical 
flames  the  correlation  would  differ  from  Eq.  (4)  be¬ 
cause  U  and  W  are  not  relevant  parameters;  instead, 
the  degree  of  wrinkling  and  the  propagation  speed 
of  a  spherical  flame  depend  on  a  different  parame¬ 
ter:  the  total  time  during  which  the  flame  is  exposed 
to  the  turbulence.  The  best  approach  is  to  identify 
the  category-dependent  parameters  and  then  create  a 
separate  database  and  correlation  of  data  for  each  cat¬ 
egory. 


3.  Experimental  arrangement 


2.  Categories  of  premixed  flames 

The  decision  was  made  to  study  a  large  Bun¬ 
sen  flame,  which  has  been  labeled  an  “envelope 
flame”  by  Cheng  and  Shepherd  [16]  because  all  of 
the  reactants  must  pass  through  the  flame.  Other 
categories  [16]  are  “oblique  (rod- stabilized)  flames” 
which  have  been  documented  by  Knaus  and  Gouldin 
[17],  “unattached”  (low- swirl  flames  and  counterflow 
flames  [16]),  spherical  expanding  flames,  and  shear- 
dominated  (jetlike)  premixed  flames  [18].  For  each 
category  the  wrinkling  process  and  the  boundary  con¬ 
ditions  are  different.  Thus,  the  formula  for  the  burning 
velocity  for  one  category  is  expected  to  differ  from 
that  for  another  category.  The  formula  should  include 
a  few  category- specific  parameters  such  as  burner 
width  (IT)  and  mean  velocity  (U).  Experimental  evi¬ 
dence  indicates  that  there  is  no  turbulent  flame  that  is 


The  relatively  large  two-dimensional  slot  Bunsen 
burner  shown  in  Fig.  1  was  operated  under  the  con¬ 
ditions  listed  in  Table  1.  One  advantage  of  Bunsen 
geometry  is  that  all  of  the  reactants  pass  through  the 
flame  so  the  mass  flow  rate  of  the  reactants  that  are 
consumed  is  known.  The  large  burner  width  (IT)  of 
25.4  mm  creates  a  tall  flame  with  a  height  exceed¬ 
ing  100  mm;  this  allows  the  flame  to  become  highly 
wrinkled,  so  it  provides  a  good  test  of  the  ability  of 
a  model  to  predict  the  correct  degree  of  wrinkling. 

The  rectangular  central  burner  in  Fig.  1  contains 
a  turbulence-producing  grid,  and  this  central  burner 
is  surrounded  by  two  outer  burners;  all  three  burn¬ 
ers  have  the  exit  dimensions  W  =  25.4  mm  and  L  = 
50.8  mm.  Each  outer  burner  contains  a  short  grid- 
stabilized  flames  to  provide  a  stream  of  products  that 
have  the  same  axial  velocity  as  the  products  produced 
by  the  central  flame.  This  eliminates  any  interaction 
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(a) 

aX 

L 


Fig.  1.  The  2D  slot  Bunsen  flame  (a)  and  burner  (b).  Width 
W  is  25.4  mm;  length  L  is  50.8  mm.  Two  outer  burners  elim¬ 
inate  shear  layers.  P  is  the  perimeter  of  the  smooth  curve 
from  A  to  B  to  C,  obtained  by  overlapping  70  images.  Stoi¬ 
chiometric  methane-air  was  used  for  all  conditions,  (c)  Con¬ 
trol  volume  considered. 

between  the  central  flame  and  shear  layers  that  are 
associated  with  the  interface  between  products  and 
the  surrounding  room  air.  Therefore  all  of  the  flame 
wrinkling  is  due  to  the  grid-generated  turbulence.  The 
reactants  were  stoichiometric  methane  and  air,  me¬ 
tered  by  mass  flow  controllers,  for  all  cases  and  for 
all  three  burners.  The  pressure  was  1  atm,  the  reac¬ 
tants  entered  at  296  K,  Sp 0  is  38  cm/ s,  the  unstretched 
laminar  flame  thickness  <$lo  (—  7.4g'lo/Slo)  [19]  was 
0.35  mm,  and  the  Markstein  number  was  0.95  [20]. 


The  turbulence  intensity  was  varied  by  replacing  the 
central  burner  with  an  identical  burner  that  contained 
a  different  grid  or  a  different  layer  of  metal  spheres. 
Table  1  shows  the  values  of  rms  velocity  fluctuations 
u'  and  mean  velocity  £/cl,0»  which  were  measured 
with  a  laser  velocimeter  on  the  burner  centerline  just 
above  the  burner  exit.  The  grid  that  imparted  the 
lowest  turbulence  level  was  denoted  the  5%  burner; 
u'/Uql,0  varied  from  4  to  7%  for  different  values 
of  U,  but  the  average  value  was  5%.  The  other  two 
burners  were  denoted  the  10  and  20%  burners,  respec¬ 
tively.  The  spatially  averaged  burner  exit  velocity  (U) 
is  defined  as  the  measured  mass  flow  rates  of  reactants 
provided  to  the  central  burner  divided  by  the  burner 
exit  area  and  reactant  density.  Values  of  U  were  set  at 
3,  5,  8,  and  12  m/s  for  cases  3,  5,  8,  and  12,  respec¬ 
tively.  The  integral  scale  t  at  the  burner  exit  was  de¬ 
termined  in  the  standard  way  from  PIV  data  by  com¬ 
puting  the  spatial  correlation  u'(x)uf(x  +  Ax)/u'2 
and  integrating  over  all  Ax  where  x  is  the  vertical 
streamwidth  coordinate.  Integral  scales  varied  from 
3.1  to  12.8  mm  and  are  listed  in  Table  1.  For  a  typical 
condition  (case  8c),  the  turbulence  Reynolds  number 
u'l/v  was  990,  and  the  Karlovitz  number  based  on  the 
integral  scale  (ur  / 1)  /  (S^/ oiq)  was  0.020.  Bray  and 
Cant  [15]  showed  that  Iq  in  Eq.  (1)  is  1-0.28  KaMa\ 
thus  Iq  has  a  value  of  0.994. 

To  determine  if  any  of  the  reactants  can  avoid 
passing  through  the  flame,  which  is  undesirable,  oil 
drops  were  added  to  the  reactants  and  were  illumi¬ 
nated  by  a  laser  light  sheet  and  no  oil  drops  were 
observed  to  pass  out  the  end  planes  (the  plane  in 
which  points  A,  B,  and  C  are  located  in  Fig.  1).  Also, 
no  oil  drops  were  observed  anywhere  in  the  product 
gases.  Thus  it  is  concluded  that  essentially  none  of  the 
reactants  can  bypass  the  flame  sheet.  It  also  was  de¬ 
termined  that  the  mean  flow  field  is  two-dimensional 
and  does  not  vary  in  the  z  direction  (shown  in  Fig.  1) 
over  the  region  z  =  ±12.5  mm,  as  shown  in  Appen¬ 
dix  A.  There  were  some  edge  effects  of  the  end  planes 
(z  =  +25  mm  and  z  =  —25  mm)  but  all  measure¬ 
ments  were  made  in  the  x-y  plane  (z  =  0),  which  is 
far  from  the  edges,  so  edge  effects  are  not  significant. 

The  simultaneous  CH  PLIF-PIV  system  is  de¬ 
scribed  in  Appendix  A  and  in  Refs.  [21-23],  however 
many  updates  to  the  PIV  system  have  been  made;  it  is 
located  at  Wright-Patterson  AFB.  A  35  x  35-mm  re¬ 
gion  was  imaged  at  three  different  heights  above  the 
burner.  The  spatial  resolution  of  the  CH  PLIF  system 
was  0.2  mm,  which  equals  the  PLIF  laser  sheet  thick¬ 
ness  and  the  size  of  the  3  x  3 -pixel  region  in  the  field 
of  view  over  which  the  PLIF  signal  was  averaged.  The 
spatial  resolution  of  the  PIV  system  was  set  by  the 
0.27  mm  interrogating  box  size  (16  x  16  pixels)  and 
the  use  of  50%  overlap  during  PIV  data  reduction,  and 
the  size  of  the  PIV  interrogation  box  used. 
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Table  1 

Properties  of  the  12  turbulent  premixed  flame  casesa 


Case 

U 

(m/s) 

Ucl,0 

(m/s) 

u' 

uC 

SU) 

Integral  scale,  i 
(mm) 

Taylor  scale,  Xj 
(mm) 

Kolmogorov  scale,  A,r 
(mm) 

UCL,0 

3a 

3.0 

2.7 

0.04 

0.28 

- 

- 

- 

3b 

3.0 

2.7 

0.10 

0.69 

5.2 

0.5 

0.2 

3c 

3.0 

3.4 

0.24 

2.14 

3.1 

0.2 

0.07 

5a 

5.0 

4.7 

0.05 

0.68 

- 

- 

- 

5b 

5.0 

4.4 

0.10 

1.12 

- 

- 

- 

5c 

5.0 

5.9 

0.23 

3.57 

- 

- 

- 

8a 

8.0 

7.3 

0.05 

1.00 

- 

- 

8b 

8.0 

6.5 

0.11 

1.88 

12.8 

0.5 

0.1 

8c 

8.0 

8.5 

0.20 

4.47 

8.7 

0.3 

0.05 

12a 

12.0 

11.0 

0.07 

2.17 

- 

- 

- 

12b 

12.0 

11.4 

0.16 

4.80 

- 

- 

- 

12c 

12.0 

13.0 

0.25 

8.55 

- 

- 

- 

a  U  is  the  mean  velocity  that  is  spatially  averaged  over  the  burner  exit;  it  is  the  mass  flow  rate  of  reactants  divided  by  the 
burner  exit  area  and  the  reactant  density.  Ucl  0  ar*d  u'  were  measured  at  the  origin  using  laser  velocimetry.  Spo  is  0.38  m/s, 
($L0  is  0.35  mm,  burner  width  (TT)  is  25.4  mm,  and  Ma  =  0.95. 


4.  Definition  of  the  turbulent  burning  velocity: 
the  global  consumption  speed  (Sj) 

There  is  no  single  accepted  definition  of  turbu¬ 
lent  burning  velocity;  instead,  four  different  velocities 
have  been  used  in  previous  studies.  These  are  the  lo¬ 
cal  and  global  displacement  speeds  (S&  and  S^)  and 
the  local  and  global  consumption  speeds  (Sj  and  Sj). 
Each  represents  an  acceptable  way  to  characterize 
turbulent  burning  velocity,  providing  that  one  consis¬ 
tently  uses  the  same  definition  when  comparisons  are 
made.  Local  displacement  speed  S&  is  the  relative  ve¬ 
locity  between  some  contour  (such  as  c  =  0.5)  and 
the  gas  velocity  at  some  defined  distance  ahead  of 
the  flame;  it  is  difficult  to  measure  in  Bunsen  and 
V-flames  because  these  contours  are  oblique  to  the 
incoming  flow  [24].  Measurements  of  the  local  con¬ 
sumption  speed  [25,26]  require  the  measurement  of 
mass  fluxes,  including  the  correlation  p'uf,  on  all 
boundaries  of  a  small  control  volume. 

In  the  present  study,  it  was  decided  to  measure  the 
global  consumption  speed  (Sj).  If  we  consider  any 
wrinkled  turbulent  wave,  the  speed  of  the  entire  wave 
is  the  speed  of  the  smooth  surface  that  represents 
the  time-averaged  position  of  the  wave.  This  time- 
averaged  wave  position  is  identified  by  the  smooth 
curve  A-B-C  in  Fig.  lb.  The  frontal  area  of  the  time- 
averaged  flame  position  ( Ap)  is  defined  as  the  perime¬ 
ter  (P)  of  the  smooth  curve  A-B-C  multiplied  by  the 
length  L  shown  in  Fig.  1 ;  thus, 

A¥  =  PL.  (5a) 

Now  consider  the  control  volume  shown  in  Fig.  lc, 
which  encloses  the  entire  flame.  The  mass  per  second 
of  reactants  that  enter  through  the  upstream  bound¬ 
ary  of  the  control  volume  must  equal  the  known  mass 


flow  rate  provided  to  the  burner  (wr),  so  conserva¬ 
tion  of  mass  requires  that 

=  PR^T^F-  (5b) 

The  quantity  /Or  is  the  known  density  of  the  stoichio¬ 
metric  methane-air  mixture  at  296  K,  1  atm.  It  is 
appropriate  to  use  /Or  in  Eq.  (5b)  since  it  is  the  density 
of  the  unburned  gas  through  which  the  wave  is  prop¬ 
agating.  Sj  represents  the  global  consumption  speed 
because  it  is  proportional  to  the  total  mass  per  second 
of  reactants  consumed.  Eq.  (5b)  is  rearranged  to  yield 

ST  =  wr/ (prAf).  (5c) 

To  determine  P,  70  images  of  the  Mie  scattering  of 
oil  droplets  (similar  to  that  shown  in  Fig.  la)  were  ob¬ 
tained  for  each  flow  condition.  Images  were  normal¬ 
ized  to  account  for  laser  intensity  variations  and  were 
binarized;  pixels  having  intensities  exceeding  half  the 
maximum  were  set  to  a  value  of  unity  and  all  other 
pixel  values  were  set  to  zero.  The  cross-sectional  area 
of  the  reactants  in  the  laser  sheet,  which  is  the  area 
of  the  white  region  in  Fig.  la,  was  measured.  The 
shape  of  the  curve  A-B-C  was  determined  by  super¬ 
imposing  the  images  and  then  fitting  a  parabola  to  this 
time-averaged  flame  position.  It  was  required  that  the 
cross-sectional  area  under  the  resulting  curve  A-B-C 
equal  the  average  cross-sectional  area  of  the  reactants 
in  the  70  images.  This  area-matching  requirement  cor¬ 
rectly  accounts  for  large  displacements  of  the  flame 
to  the  right  or  to  the  left  of  centerline,  and  for  large 
distortions  of  the  flame  shape  as  islands  are  formed. 
Isolated  pockets  of  reactants  are  seen  at  the  top  of 
Fig.  la;  these  pockets  were  included  when  determin¬ 
ing  the  white  area  in  Fig.  la. 

Use  of  the  global  turbulent  burning  velocity  (Sj) 
offers  several  advantages  when  the  desire  is  to  as- 
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Fig.  2.  Measured  turbulent  burning  velocity,  defined  as 
global  consumption  speed  Sj.  Conditions  are  listed  in  Ta¬ 
ble  1.  For  each  curve,  mean  velocity  U  was  held  constant. 
Stoichiometric  methane-air. 


sess  the  ability  of  a  numerical  simulation  to  pre¬ 
dict  the  correct  total  consumption  rate  of  reactants. 
Ambiguities  and  errors  are  eliminated  because  it  is 
not  necessary  to  measure  the  local  velocity  compo¬ 
nents,  the  location  of  the  leading  edge,  or  the  nor¬ 
mal  to  the  flame  brush,  which  are  required  for  lo¬ 
cal  displacement  speed  measurements  [24].  Shepherd 
and  Cheng  [25]  and  Ghenai  et  al.  [26]  argued  that 
for  oblique  flames,  the  consumption  speed  is  a  bet¬ 
ter  measure  of  turbulent  burning  velocity  than  is  the 
local  displacement  speed.  Displacement  speed  has 
been  shown  to  be  three  times  larger  than  consump¬ 
tion  speed  [26]  and  is  not  a  good  measure  of  the  local 
ability  of  the  flame  to  convert  reactants  to  products 
when  the  flame  is  oblique  to  the  oncoming  flow. 


5.  Results 

5.1.  Global  consumption  speeds  Sj 

Fig.  2  is  a  plot  of  the  measured  values  of  turbulent 
burning  velocity,  which  is  defined  as  the  global  con¬ 
sumption  speed  Sj  given  by  Eq.  (5b).  Mean  velocity 
and  turbulence  levels  m'/Slo  were  varied  up  to  the  rel¬ 
atively  large  values  of  12  m/s  and  8.55,  respectively. 
Two  conclusions  that  can  be  drawn  from  Fig.  2  are 
that  the  mean  velocity  U  is  an  important  parameter  for 
this  category  of  flames,  and  that  if  U  is  held  constant 
while  u'  is  varied  independently,  the  burning  velocity 
curves  display  a  nonlinear  “bending.”  U  is  important 
because  it  affects  the  time  that  it  takes  eddies  to  cross 
the  thin  wrinkled  reaction  layer.  This  eddy  residence 


Fig.  3.  Measured  turbulent  burning  velocity  replotted  for  a 
constant  value  of  u'/U  =  0.2,  so  that  u'  and  U  are  not  in¬ 
dependently  varied.  This  method  has  been  used  previously 
and  it  incorrectly  hides  the  true  nonlinear  behavior  seen  in 
Fig.  2. 

time  is  defined  as 

fi-es  —  (U  sin  ft) ,  (6) 

where  i  is  the  integral  scale  and  ft  is  the  angle  be¬ 
tween  the  x  axis  and  the  c  =  0.5  contour,  so  U  sin  /3 
is  the  component  of  the  upstream  velocity  that  is  nor¬ 
mal  to  the  c  =  0.5  contour.  In  the  present  study,  U  was 
held  constant  while  u '  was  increased  by  using  a  differ¬ 
ent  grid,  to  increase  the  rotational  velocity  of  the  ed¬ 
dies  while  keeping  the  residence  time  approximately 
constant.  This  promotes  the  highly  wrinkled,  corru¬ 
gated  flame  conditions  that  lead  to  flamelet  merging 
and  the  “bending”  phenomenon  that  has  been  reported 
in  Refs.  [27-32]. 

It  is  reasonable  to  expect  that  “bending”  will  oc¬ 
cur,  because  as  the  turbulence  level  is  increased,  the 
flame  surface  area  cannot  grow  indefinitely.  Flamelets 
eventually  merge  and  may  extinguish  due  to  strain. 
Gas  expansion  also  creates  divergent  velocity  fields 
that  prevent  segments  from  becoming  too  densely 
packed.  When  Duclos  et  al.  [5]  increased  ur /Spo 
above  5,  the  computed  curves  of  burning  velocity  be¬ 
gan  to  bend.  Similarly,  the  bending  of  the  measured 
curves  in  Fig.  2  begins  when  w'/Slo  is  2-4.  Differ¬ 
ent  explanations  for  the  bending  phenomenon  have 
been  offered  by  Peters  [33],  Bradley  [28],  and  others 
[34-36]. 
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It  is  noted  that  this  “bending”  would  not  occur  if 
precautions  were  not  taken  to  correctly  vary  u '  and  U 
independently.  This  is  shown  in  Fig.  3.  Data  in  Fig.  2 
were  replotted  for  a  fixed  value  of  u'/U  equal  to  0.2. 
The  linear  curve  that  appears  in  Fig.  3  represents  re¬ 
sults  that  would  have  been  obtained  if  a  single  grid 
were  employed  and  only  the  mass  flow  rate  was  var¬ 
ied,  which  usually  has  been  done  in  the  past.  In  such 
studies  both  the  turbulence  level  and  the  residence 
time  of  eddies  in  the  flame  brush  incorrectly  were 
varied  simultaneously.  Thus,  as  ur  and  U  increase, 
the  eddies  become  stronger,  but  their  residence  time 
is  reduced,  and  these  two  competing  processes  can 
counteract  each  other.  This  can  mask  the  “bending” 
process  and  explain  why  it  has  not  always  been  ob¬ 
served.  The  observation  that  Sj  depends  on  U  also 
helps  to  explain  the  large  scatter  observed  when  at¬ 
tempts  were  made  to  collapse  values  of  Sj  to  a  single 
curve.  In  those  studies  the  values  of  U  were  allowed 
to  vary  randomly. 

The  global  burning  velocity  measurements  in 
Fig.  2  can  be  fit  to  the  empirical  relation 
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The  values  of  constants  B\  and  #2  that  provide  the 
best  fit  of  Eq.  (7)  to  the  data  in  Fig.  2  are  2.0  x  10-3 
and  0.16,  respectively,  and  <5lo  is  0-35  mm.  The  in¬ 
tegral  scale  is  I  and  the  thermal  diffusivity  a  of  a 
stoichiometric  methane-air  mixture  is  0.20  cm2/s. 
A  physical  explanation  that  leads  to  Eq.  (7)  is  given 
by  the  scaling  analysis  below.  The  uncertainties  in  the 
data  in  Fig.  2  are  approximately  6%  and  are  discussed 
in  Appendix  A. 

A  question  that  arises  is  whether  the  nonlinear 
trend  shown  in  Fig.  2  is  due  to  changes  in  the  sur¬ 
face  area  of  the  flame  or  to  the  internal  structure  of 
the  flamelets  as  the  turbulence  intensity  increases. 
To  investigate  this  issue,  the  perimeter  of  the  wrin¬ 
kled  flame  surface  was  measured  in  70  Mie  scatter¬ 
ing  images  and  the  averaged  value  of  the  wrinkled 
perimeter  is  Pj.  Fig.  4  is  a  plot  of  (Pj/P)2,  where 
P  is  the  perimeter  of  the  smooth  curve  that  repre¬ 
sents  the  time-averaged  flame  position.  This  nondi- 
mensional  parameter  is  one  measure  of  flame  wrin¬ 
kledness  and  it  displays  the  same  trends  as  seen  in 
Fig.  2;  the  curves  bend  as  ur  is  varied,  and  the  perime¬ 
ter  increases  with  U.  This  indicates  that  changes  in 
the  surface  area  of  the  flame  do  play  a  major  role 
in  determining  the  trends  in  the  burning  velocity  re¬ 
sults. 


0  2  4 

TURBULENCE  INTENSITY  u’/SL0 

Fig.  4.  Degree  of  flame  wrinkling,  quantified  by  the  para¬ 
meter  ( Pj/P )2.  Pj  is  the  perimeter  of  the  wrinkled  flame 
surface  that  appears  in  the  laser  sheet  images,  averaged  over 
70  images.  P  is  the  perimeter  of  the  smooth  curve  A-B-C 
in  Fig.  1  that  defines  the  time-averaged  flame  position.  The 
trend  displayed  is  similar  to  that  of  the  turbulent  burning  ve¬ 
locity  in  Fig.  2. 

5.2.  Images  of  CH  reaction  layers,  quenching  rates, 
flame-eddy  interactions 

Figs.  5-7  show  some  images  of  the  CH  reaction 
layer  and  the  vorticity  field  for  relatively  large  val¬ 
ues  of  the  turbulence  level  (V/^L0  —  4.47,  case  8c) 
when  “bending”  is  observed  in  Fig.  2.  It  is  noted  that 
pockets  of  reactants  are  seen  in  all  of  these  figures, 
and  the  CH  layers  remain  relatively  thin,  which  indi¬ 
cates  that  the  conditions  are  within  the  “corrugated 
flamelet  regime.”  Fig.  5a  is  a  magnified  view  of  a 
region  in  Fig.  5b.  An  “eddy”  is  defined  as  a  region 
where  the  magnitude  of  the  instantaneous  vorticity 
exceeds  2000  s-1;  this  value  was  chosen  to  identify 
the  main  features  of  the  flow;  less  intense  vortical  re¬ 
gions  also  exist  but  they  are  not  shown.  The  red  eddies 
are  associated  with  counterclockwise  rotational  ve¬ 
locities.  To  visualize  the  gas  velocity  near  each  eddy, 
the  measured  instantaneous  velocity  vectors  are  plot¬ 
ted  in  Figs.  5-7.  The  mean  velocity  of  the  reactants  at 
the  center  of  each  image  has  been  subtracted  from  the 
velocity  at  every  location,  which  explains  why  some 
of  the  vectors  point  downward.  Note  that  the  vectors 
represent  the  instantaneous  gas  velocity  in  the  labo¬ 
ratory  frame  of  reference,  and  not  the  flame  frame 
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Fig.  5.  Images  of  the  CH  reaction  layers  (black)  and  the 
turbulent  eddies  (red  and  blue)  for  case  8c.  U  =  8  m/s. 
MV^L0  —  4.47.  Red  and  blue  regions  have  vorticity  mag¬ 
nitude  exceeding  2000  s-1.  Stoichiometric  methane-air. 
(a)  Magnified  view  of  a  15  x  22-mm  region,  reactants  are 
on  the  right,  products  are  on  the  left;  (b)  entire  field  of  view 
(35  x  35  mm);  reactants  are  in  the  middle  region.  In  the  white 
regions  the  particle  seeding  density  was  insufficient. 

of  reference,  so  they  do  not  always  point  toward  the 
flame  on  the  reactant  side  and  away  from  the  flame  on 
the  product  side,  since  the  flame  is  not  stationary.  The 
CH  layers  are  shown  as  wrinkled  black  lines;  they  are 
the  locations  where  the  CH  PLIF  signal  exceeds  half 
the  local  maximum  signal. 

The  images  indicate  that  there  is  intense  wrinkling 
of  the  CH  layers  and  that  the  eddies  are  numerous 
on  the  reactant  side,  but  are  less  numerous  on  the 
product  side.  This  observation  is  consistent  with  the 
flame- vortex  studies  of  Mueller  et  al.  [38],  who  re¬ 
ported  that  most  of  the  vorticity  is  destroyed  as  an 
eddy  crosses  a  premixed  flame,  because  gas  expan¬ 


Fig.  6.  Another  image  of  the  CH  reaction  layers  (black)  and 
the  turbulent  eddies  (red  and  blue)  for  case  8c.  U  =  8  m/s. 
MV^L0  —  4.47.  Red  and  blue  regions  have  vorticity  magni¬ 
tude  exceeding  2000  s_1 .  Stoichiometric  methane-air.  (a,  b) 
Magnified  regions  4.5  x  4.5  mm;  (c)  full  field  of  view  of 
35  x  35  mm. 

sion  spreads  the  rotational  angular  momentum  over  a 
large  region  in  space.  Figs.  5-7  also  illustrate  the  rel¬ 
ative  roles  of  flamelet  merging  and  quenching,  which 
can  cause  the  “bending”  behavior  described  above.  It 
is  seen  that  the  CH  layers  are  continuous,  unlike  the 
discontinuous  CH  layers  reported  in  the  nonpremixed 
swirl  flame  of  Ratner  et  al.  [23].  Quenching  is  defined 
to  occur  where  the  CH  signal  on  the  layer  centerline 
is  less  than  25%  of  the  maximum  CH  signal  in  the 
image.  No  regions  of  flamelet  quenching  were  ob¬ 
served  in  any  of  the  1500  CH  PLIF  images.  Therefore 
it  is  concluded  that  in  the  present  work,  local  flame 
quenching  is  not  a  possible  reason  for  the  “bending” 
phenomenon,  and  that 

<2  =  0.  (8) 

Figs.  5-7  offer  evidence  that  merging  of  flamelets 
does  occur,  but  the  rate  of  merging  cannot  be  quan¬ 
tified  with  the  present  diagnostics.  Merging  is  about 
to  occur  wherever  there  are  two  CH  reaction  layers 
that  are  propagating  toward  each  other,  with  only  a 
thin  layer  of  reactants  between  them.  This  is  seen  in 
Fig.  5b,  where  there  is  a  “neck”  region  at  the  top 
of  the  boxed  region  and  another  neck  region  on  the 
far  right  side.  A  neck  region  is  beginning  to  form  in 
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Fig.  6  at  y  =  20  mm.  Evidence  that  merging  has  just 
occurred  is  the  existence  of  the  pockets  of  reactants 
that  are  seen  in  Fig.  6c.  Pockets  can  be  caused  only 
by  flamelet  merging.  First,  long  fingers  of  the  wrin¬ 
kled  flame  surface  are  created,  as  seen  in  Figs.  5-7, 
which  look  like  those  observed  in  the  flame-vortex 
interactions  of  Refs.  [37-40].  Then  as  the  CH  layers 
propagate  toward  each  other,  they  often  form  a  neck 
that  bums  through  to  create  a  pocket.  The  merging  of 
flamelets  and  the  creation  of  pockets  can  reduce  the 
reaction  zone  surface  area.  When  the  turbulence  in¬ 
tensity  ur  / Slo  is  reduced  from  4.47  to  less  than  one 
(cases  3a,  3b,  and  5a),  very  few  peninsulas  or  pock¬ 
ets  are  observed.  This  is  evidence  that  the  merging 
rate  increases  with  turbulence  level.  The  images  of  the 
neck  regions  in  Fig.  5  are  similar  to  DNS  of  flamelet 
merging  reported  by  Chen  et  al.  [37,  Fig.  12]  for  an 
wV^LO  value  of  10. 

The  images  indicate  that  the  wrinkles  in  the  CH 
reaction  layers  are  significantly  larger  than  even  the 
largest  eddies.  The  eddies  in  Figs.  5-7  vary  from  a 
maximum  size  of  3  mm  to  a  minimum  of  0.3  mm  (the 
detecting  limit  for  the  PIY  measurement);  the  average 
size  is  about  1-2  mm.  The  average  radius  of  curvature 
of  the  CH  layer  was  measured  to  be  12.5  mm,  which 
is  more  than  six  times  the  average  eddy  size.  This  dif¬ 
ference  is  consistent  with  previous  studies  of  flame- 
vortex  interactions  [37-40].  As  an  eddy  approaches 
the  flame,  the  gas  expansion  creates  a  dilatation  ve¬ 
locity  that  greatly  exceeds  u ',  so  an  eddy  cannot  roll 
up  the  flame  and  create  a  radius  of  curvature  compa¬ 
rable  to  its  own  radius. 

Figs.  6a  and  6b  show  counterrotating  eddy  pairs 
which  exert  extensional  and  compressive  strain  on 
the  flame.  These  are  magnified  views  of  the  regions 
denoted  A  and  B  in  Fig.  6c.  To  better  visualize  the  ro¬ 
tating  velocity  pattern,  a  constant- velocity  vector  was 
subtracted  from  each  PIV  image.  Note  that  in  Fig.  6b 
the  flame  is  curved  in  the  direction  of  the  velocity 
vectors;  it  appears  that  the  flame  is  pulled  into  the 
region  between  the  two  eddies,  imparting  a  large  cur¬ 
vature  to  the  CH  layer.  In  theory,  a  purely  compressive 
strain  field  should  reduce  the  flame  area,  as  negative 
strain  exerts  a  negative  stretch  rate;  however,  the  ed¬ 
dies  in  Fig.  6b  are  affecting  both  the  strain  rate  and  the 
curvature,  and  it  is  likely  that  the  flame  area  would 
increase  if  the  CH  layer  were  advected  in  the  direc¬ 
tion  of  the  velocity  vectors.  In  Fig.  6a  the  eddy  pairs 
exert  extensional  (positive)  strain,  similar  to  that  ob¬ 
served  in  a  counterflow  geometry.  Note  that  the  CH 
layer  is  thinner  between  the  two  eddies,  which  is  ex¬ 
pected  in  a  positive  strain  field.  These  flame-eddy 
interactions  in  a  fully  turbulent  flame  are  similar  to 
the  single  flame-vortex  experiments  of  Mueller  et  al. 
[38]  and  Sinibaldi  et  al.  [39]  and  the  simulations  of 
Najm  and  Wyckoff  [40],  In  Fig.  7a  a  large  eddy  exists 


Y  (mm) 


Fig.  7.  Images  of  the  CH  reaction  layers  (black)  and  the 
turbulent  eddies  (red  and  blue)  for  case  3c.  U  =  3  m/s. 
m75lo  =  2.14.  Red  and  blue  regions  have  vorticity  mag¬ 
nitude  exceeding  2000  s-1.  Stoichiometric  methane-air. 
(a)  Magnified  view  of  a  4.5  x  4.5-mm  region;  (b)  full  field 
of  view  of  a  35  x  35 -mm  region. 

within  an  elongated  flame  wrinkle.  However,  in  most 
of  the  elongated  flame  wrinkles,  no  such  large  eddy 
is  seen.  The  fact  that  the  large  wrinkles  do  not  cor¬ 
relate  with  the  presence  of  large  eddies  nearby  may 
be  due  to  the  possibility  that  a  large  eddy  previously 
created  the  wrinkle  but  now  has  disappeared.  Another 
possibility  is  that  the  eddies  trigger  hydrodynamic  or 
thermodiffusive  instabilities  in  the  flame,  which  have 
their  own  characteristic  wavelengths. 

5.3.  Stretch  rates,  stretch  efficiency  factor  /k, 
quenching  rates 

The  stretch  rate  ( K )  along  the  CH  layers  is  the 
sum  of  the  strain  rate  Ks  and  the  curvature  contribu¬ 
tion  Kc  [11]: 

K  =  KS  +  KC,  (9) 

£s  =  -n-(n-  V)VG  +  V.VG,  (10a) 

kc  =  sl/rc.  (10b) 

To  measure  the  in-plane  component  of  K ,  the  isoline 
of  the  maximum  CH  signal  was  identified  along  with 
its  normal.  This  isoline  is  located  along  the  center  of 
the  thin  black  CH  layers  seen  in  Figs.  5-7.  Then  the 
leading  edge  of  the  flame  was  defined  as  a  second  iso¬ 
line  that  is  0.4  mm  closer  to  the  reactants  than  the  first 
isoline,  and  the  in-plane  components  of  gas  velocity 
(Vg)  were  measured  on  this  second  isoline  using  al¬ 
gorithms  described  in  Ref.  [22]  and  Appendix  A.  Ex¬ 
perimental  uncertainties  are  introduced  because 
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Case  8c 


Distance  Along  Flamefront  s  (cm) 


Fig.  8.  (a)  Instantaneous  stretch  rate  (K)  measured  in  the 
tangential  direction  (s)  along  the  wrinkled  CH  reaction  layer 
for  case  8c  using  simultaneous  PIV/CH  PLIF  diagnostics. 
Mean  stretch  rate  is  366  s-1.  (b)  Instantaneous  stretch  rate 
for  case  3b;  mean  stretch  rate  is  84  s-1 . 

(1)  only  the  in-plane  components  of  velocity  and  ra¬ 
dius  of  curvature  ( Rc )  could  be  measured, 

(2)  contributions  from  scales  less  than  the  PIV  spa¬ 
tial  resolution  (0.3  mm)  were  not  resolved,  and 

(3)  the  local  propagation  speed  has  been  replaced  by 
the  unstretched  laminar  value  in  Eq.  (10b). 

Estimates  of  these  uncertainties  are  reported  in  Ap¬ 
pendix  A. 

Fig.  8a  illustrates  the  instantaneous  stretch  rates 
K  along  the  leading  edge  isoline  for  case  8c  at 
x  =  26  mm.  It  is  seen  that  stretch  rates  oscillate  be¬ 
tween  approximately  +2000  and  —1000  s-1;  the 
mean  stretch  rate  for  this  case  (8c)  is  listed  in  Ta¬ 
ble  2;  it  is  366  s-1.  Fig.  8b  is  a  plot  of  stretch  rates 
for  case  3b  at  v  =  56  mm.  K  oscillates  between 
±400  s-1  and  its  mean  value  is  84  s-1,  which  ap¬ 
pears  in  Table  2.  For  all  cases  listed  in  Table  2,  the 
mean  stretch  rate  is  positive,  as  expected.  Note  that 
the  adjacent  data  points  in  Fig.  8  are  reasonably  well 
correlated;  there  are  typically  four  to  six  data  points 
that  are  either  monotonically  increasing  or  decreas¬ 


ing  on  each  small  segment  of  the  curve.  If  random 
errors  were  excessive,  the  adjacent  data  points  would 
not  display  this  degree  of  correlation. 

The  large  maximum  instantaneous  values  of  stretch 
rate  in  Fig.  8a,  which  exceed  1000  s-1,  may  be  due 
to  several  factors.  The  turbulent  velocity  fluctuation 
( u' )  for  case  8c  was  about  2  m/s,  and  eddies  with  a 
diameter  of  2  mm  are  observed  in  Fig.  5.  The  ratio  of 
these  numbers  is  1000  s-1.  Also,  large  stretch  rates 
were  found  to  occur  where  the  CH  layer  undergoes 
a  90°  bend.  At  these  locations  the  tangential  veloc¬ 
ity  typically  was  8  m/s  where  the  flamelet  is  parallel 
to  the  mean  flow,  but  was  nearly  zero  at  a  location 
8  mm  away  where  the  flamelet  was  perpendicular  to 
the  mean  flow.  This  change  of  8  m/s  in  an  8-mm  re¬ 
gion  leads  to  a  stretch  rate  of  1000  s_1 . 

A  stretch  rate  of  1800  s-1  is  required  [41]  to  extin¬ 
guish  a  steady,  counterflow,  stoichiometric  methane- 
air  flame  at  1  atm  with  a  reactant  temperature  of 
296  K.  There  are  some  locations  in  Fig.  8a  where 
the  instantaneous  stretch  rate  exceeds  this  value,  but 
the  CH  reaction  layer  does  not  extinguish  anywhere. 
Extinction  requires  that  a  sufficient  stretch  rate  be  im¬ 
posed  for  a  sufficient  residence  time,  but  eddies  pass 
through  the  present  flame  with  large  convection  ve¬ 
locities  that  are  10  times  larger  than  those  used  in 
flame- vortex  computations  and  experiments  [12,37— 
40],  so  the  residence  time  defined  by  Eq.  (6)  is  small. 
Donbar  et  al.  [22]  found  that  nonpremixed  turbulent 
flames  also  do  not  extinguish  where  the  strain  rate 
exceeds  the  steady- state  extinction  limit  because  of 
unsteady  effects. 

The  stretch  efficiency  factor  7k  was  determined 
by  dividing  the  measured  mean  stretch  rates  in  Ta¬ 
ble  2  by  the  values  of  u'/l  listed  in  Table  1.  Fig.  9 
shows  that  /k  increases  as  the  square  of  the  inte¬ 
gral  scale,  which  is  consistent  with  the  flame-vortex 
DNS  of  Meneveau  and  Poinsot,  as  was  shown  in 
Fig.  14  of  Ref.  [12].  Thus,  larger  eddies  in  our  tur¬ 
bulent  flame  are  more  efficient  at  stretching  the  flame 
than  the  small  eddies;  this  finding  is  similar  to  that  of 
the  flame- vortex  computations  [12].  Fig.  9  also  indi¬ 
cates  that  as  the  mean  velocity  of  the  Bunsen  burner 
is  increased,  there  also  is  an  increase  in  the  stretch 
efficiency  factor.  This  can  be  explained  by  the  fact 
that  larger  mean  velocities  force  the  flame  to  become 
taller,  and  the  farther  a  flame  segment  is  located  away 
from  the  rim  (which  anchors  the  base  and  prevents 
wrinkling),  the  more  wrinkled  the  segment  becomes. 
A  wrinkled  flame  has  a  larger  probability  that  the 
segment  will  undergo  a  sharp  bend,  which  causes 
the  velocity  that  is  tangential  to  the  flame  to  change 
abruptly,  leading  to  a  large  stretch  rate.  In  addition,  for 
a  taller  flame,  the  residence  time  of  eddies  in  the  brush 
is  larger  than  the  situation  for  which  eddies  travel  nor¬ 
mal  to  the  flame  brush.  A  reasonable  curve  fit  to  the 
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Table  2 

Measured  values  of  the  mean  stretch  rate  ( K ),  mean  strain  rate  (Ks),  mean  curvature  contribution  (Kc),  the  standard  deviations, 
and  the  stretch  efficiency  function  /k 


Case 

X 

(mm) 

K 

(s-1) 

Ks 

(s-1) 

Kc 

(s-1) 

cr 

(s-1) 

(s-1) 

(s-1) 

Stretch  efficiency 
function  /  k 

3b 

21 

79 

76 

2.7 

271 

229 

146 

1.57 

3b 

56 

86 

94 

-8.4 

349 

278 

211 

1.71 

3c 

13 

184 

165 

19 

670 

581 

333 

0.70 

3c 

31 

174 

141 

33 

448 

361 

265 

0.66 

8b 

26 

223 

226 

-3.0 

1155 

1106 

333 

4.00 

8b 

61 

248 

239 

9.1 

910 

849 

329 

4.44 

8b 

96 

229 

241 

-12 

861 

816 

275 

4.10 

8c 

26 

389 

386 

3.0 

972 

909 

345 

1.99 

8c 

61 

232 

238 

-5.7 

1155 

1107 

331 

1.19 

Normalized  integral  scale  C/8L0 


Fig.  9.  Measured  stretch  efficiency  factor  (/k)  showing  that 
large  eddies  and  a  large  mean  velocity  lead  to  more  efficient 
stretching  of  the  flame  surface.  Values  of  stretch  rate  are 
listed  in  Table  2.  Laminar  flame  thickness  ($lo  is  0-35  mm. 


data  in  Fig.  9  is 


rK  =  c  1 


(11) 


The  constant  C\  is  determined  to  be  6.6  x  10-6, 
based  on  a  laminar  flame  thickness  <$lo  (=  7.4^0/ 
^Lo)  equal  to  0.35  mm.  Fig.  10  shows  that  the  proba¬ 
bility  density  functions  (pdf)  of  the  measured  strain 
rates  Ks  and  flame  curvature  are  nearly  Gaussian 
functions.  Values  of  the  means  and  standard  devia¬ 
tions  of  the  pdf  values  are  listed  in  Table  2. 

5.4.  A  scaling  relation  that  explains  the  burning 
velocity  trends 

A  scaling  analysis  was  conducted  to  better  un¬ 
derstand  the  physical  reasons  why  burning  velocities 
in  Fig.  2  depend  on  mean  velocity  U  and  burner 
width  W ,  and  why  the  curves  display  “bending.”  It  is 
now  shown  that  Eq.  (7),  which  is  an  empirical  curve 
fit  to  the  data,  can  be  derived  from  the  conservation 
equation  for  flame  surface  density.  Consider  a  turbu¬ 
lent  Bunsen  flame  that  is  fairly  tall  so  that  the  contour 
of  mean  reactedness  c  =  0.5  is  nearly  aligned  with 
the  mean  flow  velocity  U.  This  will  be  denoted  the 
tall  flame  assumption.  The  x  direction  is  nearly  par¬ 
allel  to  the  flame  brush  and  the  coordinate  normal  to 
the  brush  (ij)  can  be  approximated  by  y.  The  Z  equa¬ 
tion  (Eq.  (2))  is  multiplied  by  dy  and  integrated  from 
y  =  0  to  -boo: 

OO  00 

r  -dz  r  -dz 

)u52dy  +  Jv 
0  0 
00  _  00 

r  d2z  r  - 

=  J  vT-^2 ~dy  +  J  ( K-m)Zdy .  (12) 

0  0 

In  Eq.  (12)  the  quenching  term  Q  has  been  ne¬ 
glected  because  no  quenching  is  observed,  as  stated 
by  Eq.  (8).  A  new  merging  quantity  m  has  been  de¬ 
fined  to  be  M/Z.  Because  flame  brush  is  assumed 
to  be  nearly  vertical,  the  second  term  in  Eq.  (12)  can 
be  neglected,  as  V  is  much  less  than  U.  The  third 
term  in  Eq.  (12)  is  zero,  because,  after  integration,  it 
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Fig.  10.  Probability  density  functions  of  (a)  the  strain 
rates  and  (b)  flame  curvature  for  case  8c.  U  =  8  m/s. 
w'/Slo  —  4.47.  Mean  values  and  standard  deviations  are 
given  in  Table  2. 

is  the  difference  between  dU/dy  evaluated  at  y  =  0 
and  infinity;  both  of  these  values  are  zero.  For  sim¬ 
plicity,  it  is  assumed  that  the  following  quantities  do 
not  vary  in  space:  K,  m,  the  turbulent  diffusivity  vj, 
U,  u\  and  the  integral  scale.  These  assumptions  are 
not  realistic  but  the  goal  is  to  integrate  the  equations 
under  ideal  conditions  to  show  why  certain  govern¬ 
ing  parameters  appear.  The  quantity  /  U  dy  equals 
approximately  (St/^LO)  fg-1  because  of  Eq.  (1),  and 
70  is  approximately  equal  to  unity  for  stoichiometric 
methane-air  reactants.  Therefore  Eq.  (12)  can  be  re¬ 
arranged  to  yield 

Ud(ST/SL0)/dx  =  (K-m)(ST/SL0),  (13a) 

which  is 

(St/Sid)-1  </(St/Slo)  =  [(K  -  m)/U]  dx.  (13b) 
Integration  yields 

ST/SL0  =  txp[(K  -m)x/U].  (13c) 


The  right-hand  side  of  Eq.  (13c)  is  expanded  in  a 
series  and  becomes  (1  +  [(K  —  m)x/U])  plus  higher- 
order  terms  which  are  neglected.  To  justify  neglecting 
the  higher-order  terms,  note  that  values  in  Tables  1 
and  2  indicate  that  K  =  223  s-1  when  0  =  8  m/s, 
x  =  26  mm;  the  second-order  term  in  the  expansion 
can  be  shown  to  be  less  than  25%  of  the  first-order 
term.  The  global  burning  velocity  (Sj)  is  defined  as 
the  average  value  of  the  local  burning  velocity  (St), 
so  that 

H 

St/Slo  =  H-1  j  ( ST/SL0)dx ,  (14) 

0 

where  H  is  the  average  height  of  the  Bunsen  flame. 
Combining  Eqs.  (13c)  and  (14)  yields 

St/Slo  -  1  =  (1/2 )(K  -  m)(H/U).  (15) 

The  physical  meaning  of  Eq.  (15)  is  clear;  the 
global  propagation  speed  depends  on  the  product  of 
the  effective  stretch  rate  (K  —  m)  and  the  geometric 
time  scale  ( H/U ),  which  represents  the  time  for  an 
eddy  to  travel  along  the  flame  from  the  base  to  the  tip. 
That  is,  wrinkling  is  due  to  a  sufficiently  large  rate 
of  stretching  of  the  surface  that  is  applied  for  a  suf¬ 
ficiently  long  time.  The  height  of  any  Bunsen  flame 
(. H )  can  be  approximated  by  the  geometric  relation 
(W / 2) (tan /?)“!,  where  is  half  of  the  included  an¬ 
gle  at  the  tip.  From  geometry,  Sj  is  U  sin/3,  and  for 
small  angles  sin  ft  =  tan  so 

H  =  (W/2)(U/ST).  (16) 

Now  Eq.  (16)  is  inserted  into  Eq.  (15)  to  yield 

5t/Slo  =  1  +  0.5(W7SL0)1/2(^  -  m)l/2-  (17) 

Conditions  are  assumed  to  be  consistent  with  the  tall 
flame  approximation  discussed  above,  so  the  ratio 
5t/5lo  is  much  greater  than  unity.  The  stretch  effi¬ 
ciency  function  (/k)  and  the  merging  efficiency  func¬ 
tion  (/m)  are  defined  in  the  same  manner  as  is  done 


in  Ref.  [12]: 

K  =  (u'/e)rK,  (18a) 

m  =  (18b) 

Combining  Eqs.  (17)  and  (18)  yields 

St/SLo  =  1  +  0.5(W/£)1/2(MysL0)1/2 

X  (rK  -  rM)1/2.  (19) 

As  stated  in  the  previous  section,  the  measured  values 
of  7k  in  Fig.  10a  and  10b  are  best  fit  by 

rK  =  Cl  (£/SL0)2(ty/SL0)2.  (20a) 
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The  merging  rate  m  was  not  measured  so  the  quan¬ 
tity  /m  in  Eq.  (18b)  can  only  be  estimated  using  the 
relation 


rM  =  C2(u'/Sw)(t/8w)2(U/Sw)2.  (20b) 


This  estimate  is  based  on  the  observation  that  flamelet 
merging  increases  when  the  stretch  rate  increases 
(as  the  flame  becomes  more  wrinkled),  so  it  is  log¬ 
ical  that  rM  is  proportional  to  /k-  It  also  is  noted 
that  Eq.  (20b)  can  be  justified  entirely  on  empirical 
grounds.  This  analysis  is  not  a  prediction  that  is  based 
on  first  principles;  instead,  Eq.  (20b)  seems  to  provide 
a  good  empirical  fit  to  the  measured  burning  veloci¬ 
ties.  If  the  laminar  flame  thickness  <$lo  in  Eqs.  (20a) 
and  (20b)  is  replaced  by  its  known  value  of  ao/^LO 
and  Eq.  (19)  is  combined  with  (20a)  and  (20b),  it  fol¬ 
lows  that 


This  relation  is  the  same  as  Eq.  (7),  which  represents 
a  reasonable  fit  to  the  measurements  in  Fig.  2.  The 
values  of  constants  B\  and  #2  are  2.0  x  KT3  and 
0.16,  respectively. 

Eq.  (21)  is  useful  in  that  it  identifies  the  burner 
width  W  as  a  category-dependent  parameter.  For  the 
Bunsen  burner  (envelope  flame)  category,  a  wider 
burner  leads  to  a  taller  flame,  which  becomes  more 
wrinkled  than  a  short  flame,  and  this  increases  Sj. 
Eq.  (21)  also  explains  why  the  burning  velocity  curves 
in  Fig.  2  depend  on  mean  velocity.  Increasing  U  does 
several  things:  it  reduces  the  residence  time  (Eq.  (6)) 
if  the  flame  height  is  held  constant,  but  it  also  has  the 
competing  effect  of  making  the  flame  longer,  which 
tends  to  increase  the  residence  time.  Increasing  U 
changes  the  orientation  between  the  flame  brush  and 
the  streamlines,  which  also  affects  the  eddy  residence 
time. 

The  above  analysis  indicates  that  the  “bending”  is 
due  to  two  factors:  the  merging  of  flamelets  and  a  geo¬ 
metric  factor  associated  with  Bunsen  flames.  Merging 
of  flamelets  tends  to  cause  a  decrease  in  the  flame 
surface  area  and  Sj  as  ur  increases.  The  geometric 
factor  that  causes  bending  is  associated  with  flame 
height;  increasing  the  turbulence  level  tends  to  make 
the  flame  propagate  faster,  but  this  is  counteracted  by 
the  fact  that  a  faster  Bunsen  flame  becomes  shorter 
(Eq.  (16)),  which  reduces  the  residence  time  (H/U  in 
Eq.  (15))  during  which  stretch  is  applied. 

Eq.  (21)  contains  the  factor  U  on  the  right-hand 
side,  and  this  result  is  consistent  with  the  measured 
curves  in  Fig.  2,  which  display  increasing  height  as  U 


Fig.  11.  Flame  surface  density  (E)  measured  for  (a)  case  3b 
and  (b)  case  3c.  U  =  3  m/s.  u'/S^q  =  0.69  and  2.14,  respec¬ 
tively. 

increases.  Eq.  (21)  also  indicates  that  Sj  should  in¬ 
crease  as  the  laminar  burning  velocity  Slq  increases. 
Therefore  Sj  should  be  maximum  if  the  fuel-air  ra¬ 
tio  is  nearly  stoichiometric,  which  is  realistic.  Some 
previous  formulas  for  St  are  unrealistic  because  they 
indicate  that  Sj  becomes  independent  of  Slo  at  suffi¬ 
ciently  large  turbulence  levels;  this  implies  that  the 
propagation  speed  does  not  approach  zero  as  the 
equivalence  ratio  approaches  zero,  which  is  unreal¬ 
istic. 

5.5.  Flame  surface  density,  wrinkling,  thickness  of 
flamelets 

Profiles  of  flame  surface  density  (E)  appear  in 
Figs.  11  and  12  for  four  operating  conditions.  To  de¬ 
termine  E,  it  was  assumed  that  the  average  flame 
surface  area  per  unit  volume  equals  the  average  flame 
perimeter  per  unit  area  in  the  laser  sheet.  At  each  lo¬ 
cation  of  interest,  an  interrogation  box  of  size  2.05  x 
2.05  mm  was  chosen.  Each  image  of  the  CH  sig¬ 
nal  was  binarized;  the  values  associated  with  pixels 
within  the  CH  layer  (where  the  CH  intensity  is  more 
than  half  of  its  maximum  value)  were  set  to  unity 
(black);  all  other  pixels  were  assigned  a  value  of  zero 
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Fig.  12.  Flame  surface  density  (X1)  measured  for  (a)  case  8b 
and  (b,  c)  case  8c.  U  =  8  m/s,  w'/^LO  =  1.88  and  4.47,  re¬ 
spectively.  Lowest  figure  indicates  that  as  the  interrogation 
box  size  was  reduced  from  3.08  to  1.02  mm,  no  significant 
change  in  X  resulted. 


(white).  The  average  perimeter  of  CH  layers  in  the 
box  is  the  average  area  of  the  black  pixels  within 
the  box,  divided  by  the  average  thickness  of  the  CH 
layer  at  that  location.  Fig.  12c  shows  that  varying  the 
box  size  from  1.02  to  3.08  mm  does  not  affect  val¬ 
ues  of  X  determined  in  this  way.  A  larger  box  has 


AXIAL  DISTANCE  x  (mm) 

Fig.  13.  Measured  nondimensional  wrinkling  parameter  f2, 
defined  as  (1/2)  X  dy.  Values  of  X  used  to  determine 
Q  are  plotted  in  Fig.  12.  Flames  become  more  wrinkled  in 
the  v  direction,  but  some  bending  of  each  curve  occurs. 

a  higher  probability  of  containing  a  CH  layer,  lead¬ 
ing  to  a  larger  mean  perimeter,  but  when  this  larger 
perimeter  is  divided  by  the  larger  box  area  no  change 
in  X  results.  As  expected,  the  profiles  of  X  have  the 
shape  of  two  delta  functions  near  the  burner  exit,  and 
become  Gaussian- shaped  near  the  flame  tip.  The  pro¬ 
files  in  Figs.  11  and  12  do  display  some  asymmetry, 
which  is  due  to  flow  nonuniformities  in  the  burner  in¬ 
ternal  flow  field. 

It  is  useful  to  define  a  wrinkling  parameter  Q  to 
quantitatively  determine  if  future  DNS  results  predict 
the  same  degree  of  flame  wrinkling  as  is  found  in  ex¬ 
periments.  It  has  been  shown  [10,15]  that  the  integral 
of  X  across  the  flame  brush  is  a  good  measure  of  the 
degree  of  wrinkling  of  the  reaction  zone,  so  we  define 
a  wrinkling  parameter  (F2)\ 

+oo 

n=l-jsdy.  (22) 

—  OO 

The  factor  1  /2  occurs  because  we  are  integrating  over 
two  flame  brushes — one  on  the  right  and  one  on  the 
left  side  of  the  burner.  Fig.  13  is  a  plot  of  values  of  F2 
obtained  by  integrating  the  profiles  of  X  in  Figs.  1 1 
and  12.  As  expected,  the  flames  become  more  wrin¬ 
kled  in  the  v  direction,  and  wrinkling  increases  as  u' 
is  increased.  A  larger  mean  velocity  (case  8)  yields 
more  wrinkling  because  the  flame  is  taller  and  the 
flame-vortex  residence  time  (Eq.  (6))  is  increased. 
Some  bending  of  the  curves  in  Fig.  13  is  observed. 

Another  parameter  that  must  be  modeled  cor¬ 
rectly  is  the  thickness  of  the  chemical  reaction  zones. 
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Fig.  14.  Average  thickness  of  the  CH  reaction  layers  (^ch) 
measured  in  four  flames,  in  which  the  turbulence  intensity 
MV^L0  changes  by  a  factor  of  5,  from  0.69  (case  3b)  to 
4.47  (case  8c).  Note  that  CH  reaction  layer  thickness  is  not 
significantly  affected  by  turbulence  intensity.  •,  case  3b; 
▲,  case  3c;  O,  case  8b,  A,  case  8c. 

The  CH  reaction  zones  are  the  thin  black  regions  in 
Figs.  5-7;  they  appear  to  have  a  fairly  uniform  thick¬ 
ness  of  0.6  to  0.8  mm  and  are  not  broadened  by  more 
than  a  factor  of  2  at  any  location  except  where  merg¬ 
ing  is  occurring.  Fig.  14  is  a  plot  of  the  average  thick¬ 
ness  of  the  CH  layers  (<$chX  which  is  the  full-width 
distance  between  points  where  the  CH  signal  is  half 
of  its  maximum  intensity.  For  these  four  conditions 
the  turbulence  intensity  varies  by  a  factor  of  5  (from 
w7^L0  =  0.69  to  4.47)  yet  the  layer  thickness  <$ch 
remains  at  a  constant  value  of  approximately  0.7  mm. 
Therefore,  it  is  concluded  that  there  is  no  significant 
change  to  the  measured  thickness  of  the  CH  layers 
even  when  there  are  large  changes  in  the  turbulence 
level  and  the  spatial  location.  The  uncertainty  in  the 
values  of  CH  layer  thickness  is  0.2  mm,  which  is  the 
thickness  of  the  laser  sheet  and  the  spatial  resolution 
of  each  binned  3  x  3 -pixel  region.  The  bars  plotted  in 
Fig.  14  represent  the  precision  in  the  determination  of 
the  mean  based  on  averaging  70  images. 

Based  on  the  results  in  Figs.  5-7  and  14,  it  is  con¬ 
cluded  that  flamelets  are  observed  to  exist  for  our 
conditions  (and  not  distributed  reaction  zones),  al¬ 
though  it  cannot  be  determined  if  they  are  “laminar” 
flamelets  or  “thickened”  flamelets.  Peters  [33]  defines 
a  flamelet  as  a  reaction  layer  in  which  the  gradients 
of  scalars  (such  as  CH  concentration)  that  are  nor¬ 
mal  to  the  layer  are  much  larger  than  the  tangential 
gradients.  Analysis  of  the  CH  gradients  indicates  that 
the  thin  (0.7  mm  thick)  CH  layers  in  images  such 
as  Figs.  5-7  meet  this  criterion.  It  cannot  be  deter¬ 
mined  if  the  flamelets  are  “laminar,”  which  would 
require  that  the  molecular  diffusivity  in  the  preheat 
zone  greatly  exceeds  the  turbulent  diffusivity.  It  also 
cannot  be  determined  if  the  flamelets  are  “thickened” 
by  turbulence.  Dinkelacker  [42]  provides  images  of 


flamelets  for  which  turbulence  eddies  enter  and  mod¬ 
ify  the  preheat  zone. 

If  a  distributed  reaction  zone  existed,  it  is  expected 
that  the  thickness  of  the  CH  reaction  zone  would  in¬ 
crease  with  the  turbulent  diffusivity  ( u'l ).  However, 
the  measured  mean  values  in  Fig.  14  suggest  that 
such  an  increase  does  not  occur,  even  when  u '  is  in¬ 
creased  by  a  factor  of  5.  It  is  noted  that  the  actual 
layer  thickness  may  be  somewhat  less  than  the  val¬ 
ues  in  Fig.  14,  for  two  reasons.  The  instantaneous  CH 
layer  thickness  is,  in  many  locations,  at  the  limit  of 
the  spatial  resolution  of  0.2  mm.  This  resolution  limit 
most  strongly  biases  the  thinnest  flamefront  segments. 
Thus  the  true  thickness  of  the  layers  is  expected  to 
be  less  than  the  0.7  mm  value  shown  in  Fig.  14.  An¬ 
other  reason  the  true  thickness  is  expected  to  be  less 
than  0.7  mm  is  that  the  flame  surface  is  not  normal 
to  the  laser  sheet.  Based  on  the  data  of  Shepherd  and 
Ashurst  [43]  the  actual  3D  layer  thickness  is  the  2D 
value  (i.e.,  the  value  plotted  in  Fig.  14)  multiplied 
by  the  direction  cosine  that  is  denoted  gq  and  has  a 
typical  value  of  0.65  [43].  Thus  we  expect  that  the 
average  3D  CH  layer  thickness  is  0.5  mm,  which  is 
65%  as  large  as  the  values  in  Fig.  14.  Shepherd  and 
Ashurst  also  found  that  as  the  turbulence  level  var¬ 
ied,  the  ratio  of  3D  to  2D  layer  thickness  did  not  vary 
significantly.  It  can  be  concluded  that  the  mean  thick¬ 
ness  of  the  CH  reaction  zones  in  our  flamelets  is  less 
than  0.7  mm.  A  detailed  study  of  the  structure  of  the 
flamelets  would  require  a  spatial  resolution  that  is  sig¬ 
nificantly  less  than  our  value  of  0.2  mm. 

Previous  work  also  is  in  general  agreement  with 
the  present  findings;  Shepherd  et  al.  [44]  increased 

by  a  factor  of  three  and  observed  no  change  in 
the  thicknesses  of  thermal  layers  that  they  measured 
using  Rayleigh  scattering.  Tanahashi  et  al.  [45]  in¬ 
creased  u!  by  a  factor  of  two  in  their  direct  numerical 
simulations  and  report  no  increase  in  the  computed 
reaction  zone  thickness.  For  a  nonpremixed  jet  flame, 
Donbar  et  al.  [22]  showed  that  for  jet  Reynolds  num¬ 
bers  up  to  18,600,  the  CH  layers  remained  thin  and 
were  no  thicker  than  1  mm. 

The  thickness  of  the  flame  brush  (<5p)  is  defined 
as  the  full  width  at  half-maximum  of  the  profiles  of 
crms’  where  c  is  the  reactedness.  Values  of  8j  are 
plotted  in  Fig.  15.  The  brush  thickness  increases  in 
the  v  direction,  and  it  increases  as  the  level  of  tur¬ 
bulence  increases,  as  expected.  Reactedness  (c)  was 
determined  from  the  CH  PLIF  images  by  assigning  a 
value  of  c  equal  to  zero  to  the  reactants  and  a  c  of  one 
to  the  products.  Across  the  thin  CH  reaction  layer,  the 
value  of  c  varied  in  a  linear  manner.  The  profiles  of  c 
and  cYms  are  adequately  fit  by  an  error  function  and  a 
Gaussian  function,  respectively, 

C  =  0.5(1  +erf(0>-yo)/ST»,  (23) 
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Fig.  15.  Thickness  of  the  turbulent  brush  (<$t)  measured  for 
cases  3b,  3c,  8b,  and  8c.  <5t  is  the  FWHM  of  the  profiles  of 
the  rms  fluctuations  of  the  reactedness.  Conditions  for  each 
case  are  listed  in  Table  1 . 


crms/crms,max  —  exP(  0.5((y  X))/^t)  )>  (24) 

where  yo  is  the  lateral  coordinate  of  the  center  of  the 
flame  brush.  Therefore,  using  the  values  of  <$t  plotted 
in  Fig.  15,  one  can  use  Eqs.  (23)  and  (24)  to  repre¬ 
sent  the  profits  of  c  and  c[ms  in  the  present  experi¬ 
ment. 


6.  Conclusions 


(1)  A  comprehensive  data  set  describes  previously 
unreported  properties  of  a  turbulent  premixed 
flame  that  are  especially  useful  for  the  assess¬ 
ment  of  direct  numerical  simulations  and  models, 
including  local  stretch  rates,  a  wrinkledness  pa¬ 
rameter,  degree  of  flamelet  extinction,  reaction 
layer  thicknesses,  flame  surface  density  (27),  and 
global  consumption  speed,  which  is  one  type  of 
turbulent  burning  velocity. 

(2)  The  global  turbulent  burning  velocity  (Sj)  of 
a  Bunsen  flame  displays  nonlinear  dependence 
(“bending”)  as  the  turbulence  level  (zZ/Slo)  is 
increased  above  the  relatively  large  values  of 
4-8,  if  the  burner  width  (W)  is  made  sufficiently 
large,  and  if  u'  and  U  are  varied  independently. 
The  perimeter  of  the  wrinkled  flame  (Pp)  dis¬ 
plays  similar  bending  behavior. 

(3)  Global  turbulent  burning  velocities  were  best  fit 
to  the  empirical  relation 


Sj 

A = i + si 

1->L0 


x  ’Air  M1/2m1/2 
_Slo.  -^LO  J  L  SlO  J 

which  indicates  that  two  additional  parameters 
are  required  to  correlate  Bunsen  flame  data:  the 
mean  velocity  (U)  and  burner  width  (W),  both  of 
which  affect  the  residence  times  of  flame-eddy 
interactions.  A  scaling  analysis  was  developed  to 
explain  the  observed  dependence  on  U  and  W 
and  the  nonlinear  dependence  on  u! . 

(4)  Two  of  the  several  possible  reasons  for  the  non¬ 
linear  dependence  are:  (a)  the  observed  merging 
of  flamelets  and  (b)  a  geometric  effect  that  is 
associated  with  Bunsen  flames,  which  become 
shorter  and  have  less  distance  to  wrinkle  as 
they  propagate  faster.  Images  of  CH  layers  show 
merging  events:  thin  fingers  of  flame  where  seg¬ 
ments  are  propagating  toward  each  other.  Local 
flamelet  extinction  is  not  a  realistic  reason  for  the 
bending  in  this  study,  as  extinction  was  not  ob¬ 
served. 

(5)  In-plane  components  of  the  local  stretch  rate 
along  fully  turbulent  premixed  flames  were  mea¬ 
sured;  previously  such  measurements  were  ob¬ 
tained  only  for  laminar  or  nonpremixed  flames. 
The  stretch  efficiency  function  (7k)  was  mea¬ 
sured  for  Bunsen  flames  and  is  best  fit  by  the 
empirical  relation 


rK  =  Ci 


2 


2 


This  indicates  that  large  eddies  are  more  efficient 
at  stretching  the  flame  than  small  eddies,  as  was 
predicted  by  Meneveau  and  Poinsot.  It  also  indi¬ 
cates  that  the  mean  velocity  U  is  important  and  is 
believed  to  be  because  it  affects  the  velocity  com¬ 
ponent  that  is  tangential  to  the  wrinkled  flame 
surface. 

(6)  CH-PLIF/PIV  diagnostics  were  shown  to  be  ad¬ 
equate  to  image  eddy  structure,  eddy  sizes,  and 
certain  pairs  of  eddies  which  were  observed  to 
exert  positive  and  negative  strain. 

(7)  The  measured  average  thickness  of  the  CH  re¬ 
action  layer  flamelets  does  not  change  as  tur¬ 
bulence  intensity  is  increased  by  a  factor  of  5. 
Many  pockets  of  reactants  are  seen.  This  indi¬ 
cates  that  conditions  were  within  the  “corrugated 
flamelet”  regime.  If  distributed  reaction  zones 
exist,  then  the  thickness  of  the  CH  regions  would 
be  expected  to  increase  with  turbulent  diffusiv- 
ity,  which  is  proportional  to  u',  and  this  was  not 
observed  to  occur.  The  PLIF  spatial  resolution 
of  0.2  mm  was  less  than  the  average  CH  layer 
thickness  of  0.7  mm,  but  improved  resolution  is 
needed  to  better  define  flamelet  properties. 
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Appendix  A.  Details  of  the  diagnostics,  the  data 
reduction  procedures,  uncertainties,  and 
boundary  conditions 

CH  fluorescence  was  obtained  by  tuning  the 
390.30-nm,  15-mJ  output  of  a  Nd:YAG-pumped  dye 
laser  to  the  <2i(7.5)  transition  of  the  B2 U~-X2n 
( v '  =  0,  v"  =  0)  band.  Fluorescence  was  collected 
from  the  A-X  (1,  1),  (0,0)  and  B-X  (0,  1)  bands 
at  wavelengths  from  420  to  440  nm.  The  array  of 
the  PIMAX  camera  was  170  x  170  pixels  after  3 
by  3  pixel  binning.  A  Schott  KV-418  color  glass  fil¬ 
ter  provided  greater  than  six  orders  of  magnitude  of 
rejection  at  390  nm  while  maintaining  90%  transmis¬ 
sion  at  430  nm.  It  was  imperative  to  reject  scattering 
and  emission  from  the  0.5 -pm  alumina  PIV  particles 
using  BG-1  and  KV-418  filters  and  a  short  intensifier 
gate  width  of  80  ns.  Background  emission  was  further 
reduced  by  gating  the  intensifier  microchannel  plate. 
The  thicknesses  of  all  three  laser  sheets  were  mea¬ 
sured  by  translating  a  25-pm  slit  across  the  beam  at 
the  probe  region;  sheet  thickness  is  the  distance  be¬ 
tween  the  locations  where  the  photodiode  signal  was 
25%  of  the  peak  value. 

The  PIV  system  included  an  interline  transfer 
CCD  camera  with  2k  x  2k  resolution  (Redlake  ES- 
4.0)  that  was  fitted  with  an  interference  filter  (cen¬ 
tered  at  532  nm  and  with  a  10-nm  bandpass)  and  a 
Nikon  105-mm  macro  lens  set  to  an  f-stop  of  f/8.  The 
PIV  interrogation  box  was  16  x  16  pixels.  An  Epix, 
Inc.,  interface  card  and  controlling  software  (XCAP) 
were  used  for  capturing  images.  Typically,  4  sets  of 
50  images  were  collected  for  each  condition;  frames 
were  recorded  at  1-s  intervals,  a  limitation  imposed 
by  the  ES-4.0  camera.  Accordingly,  the  PIV  lasers 
were  Q -switched  at  1  Hz,  while  the  flashlamps  were 
triggered  at  10  Hz.  Timing  for  the  lasers  and  cameras 
was  controlled  with  a  Quantum  Composer  9318  delay 
generator,  which  allows  subharmonic  cuing  and  trig¬ 
gering.  Commercial  LaVision  PIV  software  was  ap¬ 
plied  that  employs  adaptive  interrogation-region  off¬ 
setting.  Mie  scattering  diagnostics  were  used  to  vi¬ 
sualize  the  boundary  between  product  and  reactants. 
The  532-nm  Nd:YAG  laser  light  that  was  scattered 
by  micron- sized  PD-23  oil  drops  was  recorded  with  a 
high-resolution  (3k  x  2k)  Kodak  460  digital  camera. 
The  resolution  provided  by  each  pixel  was  60  pm. 


The  uncertainty  in  the  global  consumption  speed 
(A St)  was  determined  by  taking  the  logarithmic 
derivative  of  Eq.  (5c)  to  obtain 

ASt/St  =  [(Aotr/otr)2  +  (A AF/AF)2]1/2.  (A.l) 

The  mass  flow  rate  measurement  had  an  uncertainty 
(AmR/mR)  of  3%  and  the  error  in  Ap  was  5%,  so 
Eq.  (A.  1 )  indicates  that  the  error  in  determining  con¬ 
sumption  speed  was  6%. 

The  component  of  the  vorticity  vector  that  is  nor¬ 
mal  to  the  laser  sheet  ( coz )  is  defined  as  dv/dx  — 
du/dy  and  was  measured  using  the  following  8-point 
circulation  method  [46],  which  has  second-order  ac¬ 
curacy: 

Ui-lj-l  +2u,-j_i  +U,-+1J_1 

COj  j  =  - 

l'J  8Ay 

Yi+l,y— 1  +  2vj+i,  j  +Vj+ij+i 

8Av 

ui+l,7+l  +  2uU+ 1  +uz  — 1,;+1 

8Ay 

v,--u+1+2v,--u+v, 

8A* 

where  Av  =  Ay  =  8  pixels  (0.14  mm),  with  50% 
overlapping  of  16-pixel  interrogation  regions  (IRs). 
This  method  has  less  random  error  than  a  central  dif¬ 
ference  formula,  because  it  is  equivalent  to  central  dif¬ 
ferencing  a  velocity  field  smoothed  with  a  3  x  3  ker¬ 
nel  [46].  Random  errors  in  the  vorticity  measurements 
are  due  to  (1)  the  error  in  determining  the  peak  loca¬ 
tion  of  the  PIV  correlation  function,  (2)  contributions 
from  small-scale  eddies  that  are  not  resolved,  (3)  ther- 
mophoretic  velocity  bias,  (4)  particle-image  distor¬ 
tion  from  imaging  through  the  flame,  and  (5)  perspec¬ 
tive  distortion  error.  Error  in  determining  the  peak  lo¬ 
cation  will  depend  on  the  signal-to-noise  ratio  (SNR) 
of  the  “signal”  peak  in  the  correlation  map  and  is  of¬ 
ten  the  dominant  error.  With  a  recursive  algorithm, 
such  as  the  one  employed,  IR  shifting  will  retain  the 
maximum  number  of  particle  pairs  (and  thus  the  SNR 
will  be  maximized),  and  the  effective  particle  image 
shift  should  be  ^0.5  pixel.  With  a  numerical  simu¬ 
lation,  Raffel  et  al.  [46]  show  the  dependence  of  the 
peak-finding  rms  uncertainty,  <$dis>  on  particle-image 
shift  for  an  IR  size  of  32  x  32  pixels,  mean  particle- 
image  diameter  of  2.2  pixels,  8-bits  particle-image 
digitization,  and  various  numbers  of  particle  pairs. 
Applying  the  analysis  of  Raffel  et  al.  to  the  present 
conditions,  with  an  average  of  5.2  particle  pairs  per 
IR,  the  rms  uncertainty  vs.  particle  image  shift  in¬ 
creases  from  ^dis  <  0.005  pixel  with  no  shift  to  about 
0.04  pixel  at  0.5  pixel  image  shift.  It  should  be  noted 
that  in  the  cold  reactants,  the  typical  number  of  par¬ 
ticles  was  much  greater  than  5  (and  thus  one  would 
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expect  lower  <5dis)>  whereas  in  the  hot  products,  the 
number  of  particle  pairs  was  <  5. 

Using  an  estimated  value  of  0.5  pixel  image  shift, 
results  [46]  indicate  that  <$dis  =  ±0.04  pixel;  with  a 
typical  particle  displacement  of  8  pixels,  this  yields 
a  relative  uncertainty  of  8u/u  =  <$dis/8  ~  ±0.005  in 
the  velocity.  Because  random  errors  add  in  quadra¬ 
ture,  the  random  uncertainty  in  the  vorticity  is 

8co  =  (\/8  ±  4  •  22/8)  •  <$dis/(Av  •  A t) 


the  CH  countour  peak  within  half  a  pixel  at  each  cam¬ 
era  pixel,  then  we  might  expect  the  error  in  the  flame 
surface  normal  angle  to  be  8cp  ~  ±0.1  radian. 

In  addition  to  the  errors  sources  noted  for  the  de¬ 
rived  vorticity,  those  contributing  to  error  in  the  strain 
rate  measurements  are  due  to  uncertainties  in  the  sur¬ 
face  normal  and  the  out-of-plane  velocity  gradients. 
The  random  error  contribution  to  one  of  the  velocity 
derivatives,  dv/d y,  for  example,  from  uncertainty  in 
PIV  peak  finding  is 


^0.6Sdis/(A*  •  A0, 


where  At  is  the  separation  in  time  of  the  PIV  pulses 
(A£  =  40  and  15  ps  for  cases  3  and  8,  respectively). 
Thus,  using  these  numbers,  the  uncertainties  in  vor¬ 
ticity  are  estimated  to  be  75  and  200  s-1  for  cases  3 
and  8,  respectively;  this  represents  the  “noise  floor” 
for  measurement  of  vorticity.  By  using  the  8-point  cir¬ 
culation  method,  the  spatial  resolution  of  the  vorticity 
measurement  is  32  pixels  on  a  side  (0.55  mm)  and 
~  0.3  mm  thick.  It  is  clear  from  the  above  equation 
the  amplification  effect  of  using  smaller  IRs  on  the 
derivative  quantities  (since  Av  becomes  smaller);  the 
use  of  smaller  IRs  can  also  amplify  the  uncertainty  of 
derivative  quantities  because  of  the  reduction  in  the 
correlation  SNR. 

The  total  strain  rate  on  the  flame  is  given  by 
Eq.  (10a);  the  measured  in-plane  contribution  is 

z  du  ( du  dv\ 

^s  =  (1-nd^-n,ny(-  +  -) 

+  (1_w5)^;-  (A-3) 

The  velocity  derivatives  were  determined  from  the 
second-order  accurate  differencing  relation,  which  for 
the  v  derivative  is 


dv 

dy 


1 

—  {Vi  +  l,j-Vi-lJ+Vi  +  lJ-l 

-Vj-lJ-1  +V/+1J+1  -Vj_lj+l). 

(A.4) 


The  spatial  resolution  of  the  strain  measurement  is 
the  same  as  that  of  the  vorticity  measurement:  0.55  x 
0.55  x  0.3  mm3.  In  order  to  measure  the  normal  vec¬ 
tor  (n)  and  the  radius  of  curvature  of  the  CH  reaction 
layers,  a  third-order  polynomial  was  fit  through  every 
five  points  (each  point  separated  by  0.2  mm,  the  pixel 
resolution  on  the  CH  ICCD  camera)  along  the  isoline 
of  maximum  CH  signal.  At  the  center  of  these  five 
points,  the  curvature  (1  /Rc)  is  given  by 


l/Rc  =  V  •  n  =  T  •  (dn/ds).  (A.5) 

T  is  a  unit  vector  that  is  tangential  to  the  isoline  of 
maximum  CH  signal,  5  is  the  coordinate  along  the 
CH  isoline.  If  we  simply  assume  that  we  can  locate 


<$![-  =  V6/6  •  5(jis/ (Ay  ■  At)  ss  0.4Sdjs/(A;y  ■  Af), 
dy 

equaling  ~  50  and  130  s-1  for  cases  3  and  8,  re- 
spectively.  The  sensitivity  to  uncertainty  in  the  flame- 
surface  normal,  using  the  in-plane  strain  equation,  is 


( du  dv\ 

+  2sin<pcos<p-  I  —  -  —  I. 

Considering,  for  simplicity,  a  flame  segment  that  is 
vertically  oriented  (yp  =  0),  yields 


Here,  the  random  uncertainty  has  both  a  “floor”  com¬ 
ponent  (50  and  130  s-1)  and  a  component  propor¬ 
tional  to  the  velocity  derivatives.  If,  say,  (fy  +  f§)  = 

300  and  1000  s-1,  then  8KS  =  60  and  165  s-1,  yield¬ 
ing  percentage  random  errors  of  20  and  17%,  for 
cases  3  and  8,  respectively. 

Other  contributions  (noted  above)  to  the  error 
include  those  from  the  small-scale  eddies,  ther- 
mophoretic  velocity  bias,  flame-induced  image  dis¬ 
tortion,  and  perspective  imaging  error  (where  out- 
of-plane  particle  displacements  appear  as  in-plane 
displacements).  Error  is  introduced  by  small-scale 
eddies  because  of  the  finite  resolution  of  the  PIV 
measurement:  ~  0.27  mm  for  velocities.  In  the  reac¬ 
tants,  0.27  mm  was  less  than  10  and  3%  of  the  integral 
scales  for  cases  3  and  8,  respectively.  Relative  to  the 
Taylor  and  Kolmogorov  scales,  the  PIV  resolution  is 
comparable  to  the  Taylor  scale  (Table  2)  and  three 
to  six  times  larger  than  the  Kolmogorov  scale.  Wyn- 
gaard  [47]  integrated  the  area  under  the  spectrum  of 
scalar  gradients  and  showed  that  if  the  resolution  is 
three  times  the  Kolmogorov  scale,  the  resulting  error 
in  the  gradient  is  20%.  This  is  the  estimated  error  in 
the  strain  rate  due  to  small  scale  eddies  in  the  present 
work. 

Thermophoresis  is  the  tendency  for  a  particle  to 
drift  “down”  the  temperature  gradient.  Sung  et  al.  [48] 
found  that  the  thermophoretic  drift  velocity,  equaling 
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the  difference  between  the  gas  velocity  and  the  par¬ 
ticle  velocity,  was  ~  0.1  m/s.  If  it  is  assumed  that 
this  bias  error  occurs  over  the  flamefront  thickness 
of  ~  0.7  mm,  then  the  magnitude  in  the  bias  error 
in  the  velocity  derivative  normal  to  the  flame  sur¬ 
face  can  be  reasonably  large,  ~  140  s-1.  However, 
derivatives  away  from  the  flame  surface  or  along  the 
flame  surface  (as  with  flamefront  strain)  are  mini¬ 
mally  affected  by  thermophoresis  since  the  tempera¬ 
ture  gradient  is  small.  While  deriving  a  contribution 
to  the  uncertainty  (particularly  random  uncertainty) 
from  flame-induced  image  distortion  is  difficult,  we 
can  more  easily  estimate  the  effect  of  perspective 
imaging  error.  If  we  assume  that  the  out-of-plane  ve¬ 
locity  component  is  equal  to  the  turbulence  inten¬ 
sity  and  that  w'  =  v'  (assuming  isotropic  turbulence), 
then  at  a  10-mm  radius  from  the  image  center  with 
the  lens  240  mm  from  the  object  plane  and  a  turbu¬ 
lence  intensity  of  20%,  this  will  appear  as  an  in-plane 
displacement  of  ~  (8  pixels  x  0.20)  x  (10/240)  = 
±0.07  pixel.  Of  course,  for  this  error  to  contribute 
the  random  uncertainty  in  strain  or  vorticity,  it  would 
have  to  vary  over  a  spatial  scale,  e.g.,  the  characteris¬ 
tic  eddy  size.  Only  when  the  characteristic  eddy  size 
is  less  than  1  mm  will  perspective  imaging  error  be 
a  significant  contributor  to  random  uncertainty  of  the 
derivative  quantities. 

The  out-of-plane  contribution  to  the  strain  rate  was 
not  measured;  it  is  the  difference  between  Eqs.  (10a) 
and  (A.3),  which  is 

Ks,op  =  -(l-nj)(dw/dz) 

—  nzny(dw /dy  ±  dv/dz) 

—  nznx(dw/dx  ±  du/dz).  (A.6) 

This  out-of-plane  contribution  is  argued  to  be  be¬ 
tween  25  and  33%  of  the  total  strain  rate.  Each  of  the 
velocity  derivatives  in  Eq.  (A.6)  is  estimated  to  have 
an  equal  magnitude,  on  average,  since  grid-generated 
turbulence  is  nearly  isotropic.  If  the  normal  to  the 
flame  sheet  also  is  assumed  to  be  oriented  isotropi¬ 
cally,  then  =Uy  =n^  =  1/3,  since  the  sum  of  n \ , 
Uy  and  must  add  to  unity,  and  the  correlation  terms 
such  as  nxny  would  be  zero.  Inserting  these  values 
into  Eqs.  (A.6)  and  (10a)  provides  an  estimate  of  the 
out-of-plane  strain  rate  that  is  one-third  of  the  total 
strain  rate.  A  more  probable  flame  orientation  is  not 
isotropic,  but  is  one  that  for  which  each  segment  of 
the  flame  is  more  likely  to  be  inclined  at  about  45 
degrees  to  the  flow,  as  is  the  flame  brush.  For  this  ori¬ 
entation,  n \  and  Uy  are  equal  (and  each  has  a  value 
of  1/2)  but  nz  is  nearly  zero.  The  correlation  nxny 
is  0.5  for  this  orientation.  Inserting  these  values  into 
Eqs.  (A.6)  and  (10a)  provides  an  estimate  of  the  out- 
of-plane  strain  rate  that  is  one-fourth  of  the  total  strain 
rate. 


LATERAL  DIRECTION  Z  (cm) 

Fig.  16.  Boundary  conditions  at  burner  exit  (x  =  0)  indi¬ 
cating  the  degree  of  two-dimensionality  in  the  mean  flow 
properties.  The  z  direction  is  perpendicular  to  the  laser  light 
sheet,  as  shown  in  Fig.  1.  All  data  were  obtained  at  x  =  0, 
y  —  0. 
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Fig.  17.  Axial  velocity  and  velocity  fluctuations  along  the 
burner  centerline  (y  =  0,  z  =  0)  for  four  cases. 

Error  is  introduced  because  small-scale  eddies  are 
not  resolved  by  the  spatial  resolution  of  the  PIY  sys¬ 
tem,  which  was  0.3  mm.  This  dimension  was  less 
than  10%  of  the  integral  scale  for  case  3  and  was 
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less  than  3%  of  the  integral  scale  for  case  8.  The  Tay¬ 
lor  scales  (i (Vf/v)-1/2)  and  the  Kolmogorov  scales 
(. tiu'i/v r3/4)  are  listed  in  Table  1.  The  PIV  spa- 
tial  resolution  was  about  equal  to  the  Taylor  scale  and 
was  three  to  six  times  the  Kolmogorov  scale.  Wyn- 
gaard  [47]  integrated  the  area  under  the  spectrum  of 
scalar  gradients  and  showed  that  if  the  resolution  is 
three  times  the  Kolmogorov  scale,  the  resulting  er¬ 
ror  in  the  gradient  is  20%.  This  is  the  estimated  error 
in  the  strain  rate  due  to  small  eddies  in  the  present 
work.  Another  potential  source  of  uncertainty  is  ther¬ 
mophoresis,  which  is  the  tendency  for  a  particle  to 
drift  down  a  temperature  gradient.  Sung  et  al.  [48] 
found  that  the  thermophoretic  drift  velocity  (equal¬ 
ing  the  difference  between  the  gas  velocity  and  the 
particle  velocity)  was  0.1  to  0.15  m/s.  If  it  is  as¬ 
sumed  that  this  error  occurs  over  a  distance  of  2  mm 
(which  would  require  a  temperature  gradient  change 
of  ~  2000  K  over  this  distance),  then  the  error  in  the 
velocity  derivative  would  be  50  to  75  s-1 . 

Figs.  16  and  17  quantify  the  boundary  conditions 
at  the  burner  exit  that  are  needed  to  numerically  sim¬ 
ulate  the  experiment.  Uql,0  and  u'  were  measured 
with  a  laser  velocimeter  and  normalized  by  U,  which 
is  proportional  to  the  measured  mass  flow  rate.  Fig.  16 
indicates  that  the  flame  has  2D  mean  properties  near 
z  =  0,  which  is  the  laser  sheet  location,  as  the  exit 
conditions  deviate  by  no  more  than  10%  over  a  2-cm 
distance  in  the  z  direction. 
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